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IN-SITU  MICROWAVE  MESUREMENTS  OF  LOSSY  DIELECTRICS 

R.  J.  King,  C.  D.  Kim  and  J.  B.  Beyer 
Department  of  Electrical  and  Computer  Engineering 
University  of  Wisconsin,  Madison,  WI  S3706 

INTRODUCTION 

The  prime  objective  of  this  work  has  been  to  develop  mi;1 
instrumentation  systems,  and  measurement  and  inversion  techniques 
determining  the  complex  dielectri<  constant  f  lossy  media  in-s  tn. 
commonly  used  methods  are  not  m-situ,  e.y . ,  laboratory  m-tho-ls  r--qm  • 

removing  the  sample  from  location  mj  prepjt  .ng  if  for  testing  t.  and 

fitting  it  into  a  waveguide,  cavity,  etc  .  .>iH.h  preparation  can  si g  Meat  t  y 
alter  the  electrical  properties  by  changing  the  density,  releasing  stress 
and/or  changing  the  moisture  content.  Other  in-situ  methods  rely  on  the 
reflection  of  the  wave  from  the  surface  of  the  medium.  These  methods  given 
reasonable  accuracy  for  low-loss  media,  but  the  accuracy  suffers  when  the 
medium  has  a  large  complex  dielectric  constant,  especially  when  the  loss 
tangent  is  not  small  compared  to  unity. 

Unique  features  of  this  work,  compared  to  other  studies,  are  that  tn- 
complex  propagation  constants  (y  -  a+jfs)  of  traveling  waves  are  measured  ..sing 
phase-sensitive  homodyne  instrumentation  systems.  This  phase-sens  if  v 
property  permits  measuring  the  variations  i-f  the  ph.-  ••  well  ...  •  ►.,> 

amplitude  of  the  propagating  field.  Lomodyne  .us  ar-  ' 1  ‘■•’.own  *r  •  n.- 1  • 
linearity,  extremely  good  sensitivity  ana  great  dynamic  •  .u 

The  relationship  between  the  ..implex  propagation  constant,  . ,  o  i  t  • 

complex  dielectric  constant  e  (=  e’-jc'  )  is  determi  ne.;  ny  !  ne  .net  ;  -d  i  , 

couple  the  electromagnetic  field  with  the  medium.  Ihese  coupling  netim  and 
the  solutions  of  the  corresponding  inversion  problems  t.o  find  ►.  are  the  mam 
topics  of  this  work. 

An  invasive  method  and  a  nor.  invr .,  i  ve  method  w.  •  e  •  i.  L  ' 


invasive  method,  a  small  modulated  scattering  probe  is  inserted  into  a  test 


hole  in  the  medium.  When  the  medium  is  illuminated  by  a  spherically  expanding 
wave  originating  from  a  horn  source  outside  the  medium,  the  complex 


propagation  constant,  y,  is  measured.  After  correcting  for  spherical 
expansion,  e  can  be  calculated  in  a  straightforward  manner.  Another  advantage 
of  this  methid  is  that  y(z)  can  be  measured,  and  hence  e(z)  can  be  found,  for 
materials  which  are  graded  gradually  in  the  direction  normal  to  the  surface. 
The  chief  disadvantage  of  the  method  is  that  a  hole  is  required.  The 
dielectric  constant  of  Portage  sand  with  various  moisture  contents  were 
measured. 

The  non-invasive  method  uses  an  open-walled  rectangular  waveguide  laid  on 
the  surface  to  couple  the  traveling  wave  in  the  waveguide  with  the  medium. 
The  instrumentation  and  techniques  for  measuring  y  for  the  traveling  wave  work 
quite  well,  but  the  task  of  inverting  the  data  to  find  e  is  difficult  owing  to 
the  complicated  boundary  value  problem  which  must  be  formulated  and  solved. 
An  approximate  solution  has  been  obtained  which  correctly  shows  the  general 
trends,  but  further  work  is  needed  to  attain  reasonable  accuracy.  Results  are 
given  for  several  materials  such  as  sand,  rock  and  concrete.  Besides 
improving  the  approximations  in  the  analysis,  a  particularly  promising 
improvement  in  the  apparatus  entails  filling  the  waveguide  with  a  low-loss 
dielectric  such  that  the  dielectric  constant  of  the  medium,  relative  to  that 
in  the  waveguide,  is  scaled  down.  This  has  the  additional  advantage  of 
reducing  the  size  of  applicator  waveguide. 

Further  studies  are  suggested  which  use  two-conductor  traveling  wave 
applicators  for  coupling  the  wave  to  the  medium.  The  advantage  of  these, 
compared  to  the  rectangular  waveguide,  is  that  the  operating  bandwidth  is 
substantially  increased.  Other  instrumentation  improvements  are  suggested  for 
increasing  the  sensitivity,  dynamic  range,  accuracy  and  ease  of  use. 
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I.  THE  MICROWAVE  HOMODYNE  SYS f EM 

In  this  work,  a  microwave  humodyne  system  has  been  used  to  measur- 
attenuation  and  phase  of  the  field  to  obtain  f he  complex  propagation  con  '  >  r 
Y  (=a+je),  from  which  the  complex  dielectric  constant  of  the  medium  - 

found. 

First,  the  basic  principles  of  the  instrumentation  as  shown  in  Fig.  1  are 
described  briefly.  Further  details  on  homodyne  systems  in  general  are  given 
in  the  book  by  King  [1978]. 


ATTEN.  HYBRID  TEE  MEDIUM  | 


Fig.  1.  Microwave  homodyne  system  for  measunny  ine 

complex  propagation  constant  of  a  traveling  wave. 

The  output  of  the  CW  microwave  (200  mW  at.  4.9  GHz)  source  is  split,  at 
port  1  and  fed  into  the  two  channels:  the  reference  channel  a  ni 

information  channel.  The  CW  reference  signal,  A,  of  6  -  10  dBm  is  derived  by 

\ 


reflection  from  a  short  at  port  3  as  shown  by  the  solid  signal  flow  path. 
This  level  is  sufficient  to  ensure  that  the  mixer  operates  as  a  linear 
detector.  The  short  at  port  3  can  be  moved  to  vary  the  phase  of  the  reference 
signal.  A,  although  its  position  is  generally  fixed  in  actual  operation. 

In  the  information  channel,  the  CW  signal  which  passes  from  port  2 
propagates  through  the  unknown  medium  and  is  backscattered  by  a  modulated 
scatterer  which  actually  is  a  small  dipole  (electric  or  magnetic),  switched  at 
the  audio  angular  frequency  w  (f  =  10  KHz).  The  double  sideband  amplitude 
modulated  (DSBAM)  backscattered  signal,  b,  is  fed  into  the  mixer  via  port  2 
and  port  4,  as  shown  by  the  dashed  signal  flow  path.  This  modulated  signal  at 
the  mixer  is  of  the  form 


) b J [1  +  m  cos(<jomt)]cos(wt+<t.)  (1) 

where  m  is  the  amplitude  modulation  index  and  $  is  the  phase  difference 
between  the  reference  signal  and  the  backscattered  signal. 

The  detected  homodyne  output  signal  at  wm  is  given  by 

eout  *  Khlbl,n  C0SU)  c°s(wmt)  f°r  b/A  «  1  (2) 

where  is  the  mixer  conversion  constant. 

Thus,  the  variation  of  amplitude  |b|  and  phase  $  can  be  measured  as 
functions  of  the  position  of  the  modulated  scattering  probe.  This  information 
is  then  used  to  directly  determine  the  propagation  constant  y  of  a  traveling 
wave,  owing  to  the  phase  sensitive  nature  of  homodyne  detection.  In  contrast, 
a  long-standing  indirect  technique  for  measuring  y  has  been  to  create  a 
standing  wave  in  the  field  which  is  coupled  to  the  medium  and  to  use  video 
detection  to  measure  the  standing  wave  pattern.  But  if  the  losses  of  the 
medium  are  appreciable,  the  minima  of  the  pattern  are  not  sharp  and  spaced  a 


half  wavelength  apart.  It  then  becomes  necessary  to  analyze  i  he 
graphically  or  by  computer  to  rind  a  and  a  [Gardoil,  1975].  Furtheimoi 
is  not  always  possible  to  create  a  standing  wave  in  the  medium  in-si  ui. 

This  study  uses  the  phase  sensitive  property  of  homodyne  dete 
throughout  to  measure  y  of  various  traveling  waves  which  are  coupled  t 
medium  under  investigation. 
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II.  MEASUREMENT  METHODS 

Two  different  methods,  an  invasive  method  and  a  rion-invasive  method,  were 
studied  to  determine  the  complex  dielectric  constant. 

1.  The  Invasive  Method 

In  the  invasive  method,  the  complex  propagation  constant  of  a  spherically 
expanding  wave  propagating  in  the  medium  is  measured  by  inserting  a  modulated 
probe  into  the  medium,  as  shown  in  Fig.  2. 


Fig.  2.  Measurement  set  up  for  the  invasive  method. 

The  probe  is  a  small  loop  terminated  with  two  PIN  diodes  which  are  switched  at 
10  KHz  via  resistive  leads,  as  shown  in  Fig.  3. 

The  TEM  wave  is  launched  by  a  horn  antenna,  and  propagates  perpendicular 
to  the  surface  of  the  medium.  The  magnetic  field  inside  the  medium  is  a 
spherically  expanding  wave 
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where  z  is  assumed  to  be  in  the  radiation  zone. 
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Fig.  3.  Modulated  loop  scatterer. 


Then  the  backscattered  signal,  b,  which  is  proportional  to  the  square  of 
the  incident  field  at  the  scatterer  [King,  1978]  is  of  the  fcrm 

b  .  b  e“2£  e-j2(u  (4) 

z 

Thus,  from  (2),  the  amplitude  of  the  homodyne  output  at  wm  is 

v  t  ~  Kh ! b ! m  cos(*) 

-2aZ  (5) 

=  V0  Y  ccs(2ez) 


1 
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Therefore,  to  yet  a  the  effect  of  the  i~ ^  spreading  factor  must  be 
corrected.  To  do  this  requires  determining  the  effective  origin  of  the  horn 
antenna.  This  can  be  done  by  matching  the  measured  intensity  versus  distance 

O 

curve  for  free  space  with  the  i~c  curve. 

The  factor  of  cos(20z)  gives  zero  crossings  for  every  quarter-wavelength 
in  the  medium.  Thus,  0  can  be  easily  calculated  by  measuring  the  distance 
between  zeros.  In  practice,  several  zeros  are  traversed,  and  the  results 

averaged. 

A  typical  plot  of  (5)  which  pertains  to  the  field  in  wet  Portage  sand  is 
shown  in  Fig.  4. 

The  moisture  content  of  the  uppermost  sand  is  relatively  constant  with 
depth.  However,  as  the  water  tends  to  settle  the  moisture  content  increases 

sharply  at  a  depth  of  about  8  cm,  such  that  the  space  between  the  sand 

particles  is  filled  with  water.  Thereafter,  the  moisture  content  does  not 

increase  further.  The  wave  attenuation  therefore  increases  sharply  at  the 
boundary  where  the  sand  is  saturated  with  water,  as  shown  by  the  amplitude  of 
the  last  two  peaks  in  Fig.  4. 

The  relation  between  the  dielectric  constant,  e,  and  a  and  0  is 

0-ja  =  w/ye 

=  k0  P~r  (6) 

for  non-magnetic  material  where  kg  is  the  wave  number  of  free  space. 

Thus, 


I 


-y~ 


Fig.  4.  Typical  output  ot  the  homotlyne  system 


or. 


=  (B2-a2)/kp  (8a) 

e”  =  2aB/k2  (8b) 

Several  measurement  were  made  on  Portage  sand  with  various  moisture 
contents.  The  results  are  shown  in  Fig.  5.  The  moisture  content  used  in  this 
measurement  is  defined  by 

mass  of  the  water  contained  inn„ 

ft  -  -  ■  x  1WUa> 

w  mass  of  dry  sand 


The  disadvantages  of  this  method  are: 


2.  The  Non  -  Invasi ve  Method 


This  method  uses  a  slotted  open-walled  rectangular  waveguide  laid 
directly  on  the  surface,  as  shove.  in  Tigs,  b,  7  &  8.  The  longitudinal 
propagation  constant  of  a  traveling  wave  which  is  coupled  to  the  unknown 
medium  is  measured.  Again,  a  modulated  scatterer  and  homodyne  detection  are 
used  to  probe  the  field  of  trie  traveling  wave. 


Fig.  6.  Measurement  system  for  the  non-invasive  method. 

A  TE^q  mode  propagates  toward  the  right  in  the  closed  portion  of  the 
waveguide  in  Figs.  6  and  7.  In  the  open  portion,  this  wave  is  converted  into 
a  traveling  leaky  waveguide  mode,  where  the  fields  are  coupled  to  the  unknown 
medium.  The  right  end  of  the  waveguide  in  Figs.  6  and  7  is  terminated  with  Zq 
to  minimize  the  reflection.  The  top  wail  of  the  waveguide  is  slotted 
longitudinally  and  an  electrically  modulated  slotted  line  probe  is  inserted  to 
scatter  the  incoming  traveling  wave.  In  a  reciprocal  manner,  the 
backscattered  wave  propagates  back  to  the  homodyne  instrumentation  where  it  is 


detected. 


3- 


Since  the  bottom  wall  is  removed,  the  inside  dimension  (1.872  x  <1.9% 
inches)  of  the  open-walled  waveguide  is  same  in  width  and  slightly  greater  in 
height  than  that  of  the  closed  portion  (1.872  x  0.872  inches).  The  difference 
in  height  is  just  the  wall  thickness  (0.064  inches). 

Because  the  waveguide  is  open,  the  fields  fringe  to  the  outside  as  snown 
in  Fig.  9.  Thus,  to  confine  the  fields  to  the  medium,  flanges  of  width  w  were 
added.  The  minimum  width  was  determined  experimentally  for  the  case  where  the 
medium  is  free  space  (worst  case).  These  flanges  also  have  the  advantage  cnat 
the  boundary  value  problem  reduces  from  a  three-media  problem  to  a  two-media 
problem,  and  is  therefore  more  tractable  analytically. 

A  typical  plot  of  the  detected  signal  is  shown  in  Fig.  10.  The  field 
being  scattered  propagates  as  exp(-yz),  so  phase  sensitive  detection  gives  an 


output 


=  Vn  e'2az  cos(20z) 


;.vvr 


'  •  \  \  w  — *| 


Fig.  9.  Transverse  electric  field  distribution  of 
open-walled  waveguide  with  flanges. 
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The  attenuation  constant  a  is  easily  calculated  by  measuring  the  slop 
(-2a)  of  the  envelope  of  the  peaks  from  the  graph  when  plotted  logarithmically 
(e.g.,  in  dB).  The  phase  constant  6  is  2 u/X  where  X  is  the  longitudinal 
wavelength  in  the  guide,  since  the  zero  crossing  occurs  every  quarter 
wavelength,  0  is  easily  obtained  by  measuri  3  the  distance  between  zero 
crossings.  Again,  greater  accuracy  is  obtained  by  averaging  over  several  zero 
crossings. 

The  complex  propagation  constant  for  various  materials  (rocks,  concrete, 
sand,  styrofoam,  etc.)  were  measured  and  the  results  are  shown  in  Fig.  11.  It 
is  evident  that  the  various  media  have  a  significant  effect  on  the 
attenuation,  a,  but  the  effect  on  the  phase  constant,  0,  is  small.  For 
reference,  note  that  if  the  bottom  wall  is  perfectly  conducting  (closed 
waveguide),  a=0  and  0=0^.  The  free  space  propagation  constant,  0q,  is  also 
shown  for  reference. 


Fig.  10.  Typical  homodyne  output  for  the  non-invasive  method. 


Styrofoam 


_l _ 

0.6 


-lb- 


O 

O 

o 


Mad i .j or.  Sandstone 
Ceramic  Block 
Cement  Block 


O  S.  Dakota  Granite 
O  Indiana  Limestone 


O 


q  Cartilage  Marble 


% 

_1 _ if  » 

0.7  0.8 


- lA - > 

1.0  rad/cm 

£ 


g.  11.  Measured  a  and  e  values  for  various  materials, 
using  the  instrumentation  shown  in  Figs.  6-8. 
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III.  FIELD  ANALYSIS  OF  Tilt  OPEN-WALLED  WAVEGUIDE 
1.  Mode  Cont i gurat i on 

To  find  the  dielectric  constant  e  of  the  medium  from  the  measured  y,  the 
presence  or  absence  of  particular  field  (E  and  H)  components  must  be  known. 
This  permits  describing  the  wave  with  mode  configurations  suitable  to  an 
analytical  study,  i.e.,  solving  the  boundary  value  problem.  The  coordinate 
system  is  defined  in  Fig.  9. 

First,  the  TEX  mode  was  tried.  A  "mode  filter"  in  the  form  of  an  array 

of  thin  conducting  wires  in  x-direction,  spaced  by  1  cm,  was  placed  over  the 

open  side  of  the  waveguide.  This  filter  was  intended  to  short  out  the  Ex 

component  at  the  aperture  separating  the  interior  of  the  waveguide  and  the 

exterior  medium.  As  a  result,  Ex  was  reduced  somewhat,  but  the  measured  a  and 

&  approached  those  of  the  closed  waveguide  signifying  that  the  coupling  to  the 

medium  was  actually  reduced.  This  is  undesirable  since  the  dependence  on  z  is 

reduced.  Therefore,  the  TE  mode  analysis  was  superseded  by  the  TM  mode. 

*  y 

In  general,  any  field  can  be  decomposed  into  TMy  and  TEy  modes.  The  TMy 
mode  has  an  Ey  component  while  the  TEy  mode  does  not.  The  TMy  mode  in  the 
open  guide  couples  efficiently  to  the  TEZ  modes  in  the  closed  guide  which  only 
has  a  y-component  of  the  electric  field.  Flence,  the  TMy  mode  can  be  used  to 
describe  the  field  inside  either  the  closed  or  open  guide.  It  is  unlikely 
that  the  TEy  mode  even  exists  in  the  open  guide  since  it  is  not  excited  by  the 
Ey-component  in  the  closed  guide. 

This  argument  is  confirmed  by  measurement  which  shows  large  Hz  versus  x 
and  virtually  zero  Hy  as  given  in  Fig.  12. 

The  scatterer  used  in  this  measurement  was  a  small  loop  loaded  with  two 
tiny  pill  type  PIN  diodes  electrically  modulated  and  moved  in  the  transverse 
direction  via  a  nylon  string. 


Fig.  12.  Homodyne  output  for  the  magnetic  field  inside 
the  open-walled  waveguide  at  y  -  b/2. 


Fig.  13.  Modulated  scattering  loop  probe  for  measuring  Hz  and  H 
inside  the  open  waveguide. 


We  were  unable  to  satisfactorily  measure  any  of  the  electric  field 
components  (Ex,  E^  or  Ez)  inside  the  open  waveguide  using  an  electric 
dipole.  Coupling  with  the  metal  walls  formed  loops  such  that  the  results  were 
strongly  affected  by  the  magnetic  fields  (Hx  and  Hz)  rendering  the  results 
meaningless. 

2.  Field  Analysis  for  the  TM^  Mode 

TMy  modes  inside  the  open  waveguide  can  be  represented  by  a  magnetic 
vector  potential 


and 


T  =  u 

yv 


H  =  v  x  j 


Then,  the  scalar  wavefunction  can  be  represented  as 


OO  oo  -  Y  z 

♦(x.y.z)  =  l  sin(~)  J  ^(kv)cos(k  y)e  m  dk 

m=l  0  y  J 


■  i  *'»(=?)♦.(») 

m=l 


(10) 


where 


for  z  >  0  and  0  <  y  <  b 


-Y  z 
m 


=  l  cos ( k^y )  dy 


is  the  cosine  transform  of  i^(y).  The  dispersion  relation  is 

k^  -  (^)2  -  k2  +  y2  =  0 
0  'a  '  y  m 


(11) 


“  woEo 


5s‘ 


2 

where  kQ 


(12) 
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III 


The  wavefunctions  in  (10)  have  been  chosen  to  satisfy  the  boundary  conditions 
on  the  three  metal  walls  in  the  waveguide,  i.e.,  F  =  0  at  x^O  and  a,  and  y=0. 


The  fields  inside  the  open  waveguide  are: 


E  * 


1  3 


X  Juie^  3x3y 


(13a) 


1  y  in  _ ,mnx 


Jwe 


0  m=  1 


2?  coshr>  J  Wsin  V)e 


•r  i 
m 


dk 


1  3  i> 


J  <*>eQ  3y3z 


OO  CD 

_  1 —  V  ■— i  $in(^^-)  j  j)  (k  )k  y  sin(k  y)e 

Ja)eo  m=i  2n  a  o  y  y  m  y 


(13d) 


-Y  z 
m 


dk 


H 

x  “  '  3z 


(13c) 


*  l  -2?  /  Wk»>  cosV)e 


"V 


dk 


rn=l 


h  - 
z  "  3x 


=  I  ^cos(^)  /  ^(k  )  cos(k  y)e 

m=l  0 


-Y  z 


(13d) 


dk 


The  tactic  is  to  find  the  fields  outside  the  aperture  and  equate  the 
tangential  components  at  the  aperture  to  determine  the  equivalent  y  in  terms 
of  e.  To  this  end,  apply  the  vector  compensation  theorem  [Mittra,  1961], 
modified  to  the  case  where  one  source  is  located  inside  the  perturbed  region. 
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I  •  TJ  -  "R2  .  T7‘ )dV  =  /  (IJ  x  Tl2  -  T2  x  -R')  .  <15  (14) 


where  the  subscripts  1  and  2  denote  the  fields  excited  by  the  source  1  and 
source  2,  respectively,  and  the  primed  fields  are  the  perturbed  fields  whereas 
the  unprimed  ones  are  the  unperturbed  fields.  We  take  the  unperturbed  state 
as  a  perfect  electric  conductor  over  the  entire  y=b  plane. 


Fig.  14.  Cross-section  geometry  of  open-walled  waveguide.  In  the 

unperturbed  case,  an  electric  conductor  covers  the  aperture. 


Next,  put  a  uniform  z-directed  magnetic  line  source  in  R2  in  order  to 


find  H!  . 
lz 


\  =  uz  <5{x-x'  )5(y-y' ) 
U  2  ~  0 


which  excites  TE7  fields  with  no  z-variation.  Thus  (14)  becomes 


-Zi- 


H 


-/  6(x-x')6(y-y')H'lz(x,y,z)dx  dy  dz  =  /  (TJ  x  T?2  -  x  T(j)  •  d5  (15) 

Rj  S 

As  R  ♦  «  the  RHS  of  (15)  becomes  the  integration  over  entire  y=b  plane. 

OO  OO  CD 

/  H ^ 2 ( x '  ,y 1  ,z)dz  =  -/  {  /  (Tj  x  TT2  -  I2  x  TT* )  .•  ndxjdz  (16) 

For  the  unperturbed  state  TT  x  =  0  on  y=b  plane  (perfect  conducting  plane), 
and  for  the  perturbed  state  n  x  =  0  on  y=b  except  over  the  aperture 
(0  <  x  <  a  and  z  >  0).  Therefore,  (16)  reduces  to 


/  H'z(x',y\z)dz  -  -/  /  (l[  -  H2)  •  u  dx  dz 


■  i  /  <e;,h2z  -  e;zh2xW  «« d2 


The  field  excited  by  the  z-directed  magnetic  line  source  is  a  TEZ  mode.  Thus, 

H2x  "  am 

Since  there  is  no  z-variation  in  <j>e,  =  0. 

Moverover,  H|  for  z  <  0  is  the  magnetic  field  of  a  diffracted  wave.  Thus 
for  y'=b,  Hjz(x',b,z)  is  small  and  decreases  rapidly  for  z  <  0.  Hence,  as  a 
first  approximation,  neglect  its  contribution 


AH  =  f  H|z(x',b,z)dz 


Then  (17)  becomes 


J  H'  (x'  ,b,z)dz  =■  J  /  El  (x,b,z)H  (x,b)dx  dz  (19) 

0  0  0  1 

Both  Hjz  and  E^z  have  the  same  z-variation,  but  h2z  has  no  z-variation. 
Therefore, 


i 
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a 

Hjz(x',b,z)  =/  E^x(x,b,z)H2z(x,b)dx  (20) 

The  H2z(x,b)  of  the  source  at  (x',b)  can  be  found  by  using  the  potential 
function 

^  =  "jy  H^(kjx-x'  | ) 
and 


E  =  -V  X  U  / 

zr 


For  the  TEZ  field  of  the  line  source 


u  1  ,3  ,2.  e 

2z  j  oiu  '&z2 


-  -  2W  hS2,(Mx-x'I) 


(21) 


Then,  (20)  becomes 

H{z(x',b,z)  =  -  -g—  /  E^x(x,b,z)H^(k|x-x'  |)dx  (22) 

The  surface  impedance  as  seen  looking  from  the  waveguide  into  the  medium  is 


2  ..ik 

s  Hb 


E! 


y=b 


_lz 


y-b 


(23) 


Put  this  relation  into  (16)  and  get 


k  Z  a  /p\ 

Hjz(x'  ,b,z)  =  /  H'iz(x»'b >z )Hg  ( k | x -x *  | )dx 


(24) 


From  (13d)  H|z(x,y,z)  inside  the  waveguide  is 


Hiz(x.y.z)  -  l  ^rcos(~-)/  yk  )  cos(k  v)e 


"Y  z 
m 


dk 


m=l 


r  mit  .mux,  ,  ,, 
=  J  —  cos ( — - — )  *  (y) 
a  a  m 

m=l 


(25) 


«  t-i  8 
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M2 


(x'b)  (x,b) 


1 _ 


->  X 


Fig.  15.  Geometry  of  unperturbed  problem  for  determining  Hgj • 


Thus,  equating  the  tangential  magnetic  fields  (24)  and  (25)  at  the  aperture, 


=  cos(=il)  ♦„(»)  =4^  l  \  ^  cos(2Ii)*m(6)H(2)(Mx-x'|)<lx 
...  1  u  m=l 


I  r*»<*> 


m=I  "  0 

Equating  term-by-term  for  each  m, 


S  /  cos(— •  )ni2^  (k  I  x-x '  !)dx 


2  top 


(26) 


cos(~~)  =  /  cos(^)Hj2^(kjx-x‘ |)dx 


(27) 
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Now,  E^z  is  also  needed  in  order  to  find  the  complete  tangential  fields 
at  the  aperture.  This  makes  the  solution  unique.  The  process  can  be  repeated 
for  the  unperturbed  problem  where  source  2  is  a  uniform  electric  line  current 
and  a  perfect  magnetic  conducting  strip  covers  0  <  x  <  a  and  -«  <  z  <  +». 
Then  remove  the  perfect  magnetic  conducting  strip  to  form  the  perturbed 
problem. 

Let  72  =  uz6(x-x‘)6(y-yl) 

H2  =  ° 

which  excites  a  TMZ  field  with  no  z-variation. 


Fig.  16.  Geometry  of  the  unperturbed  problem.  A  magnetic 
conducting  strip  extends  -®  <  z  <  +®. 

Then,  (14)  becomes 

/  E|z(x' ,y ' ,z)dz  =  -/  {  /  (I2  *  TTj)y-b  •  uydx}dz  (28) 

since  H2  has  no  tangential  component  on  the  magnetic  conductor  which 
covers  the  aperture  and  E'^  has  no  tangential  component  on  the  perfect  electric 
conductors  at  the  sides  of  the  aperture. 
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Also,  for  the  TMZ  mode 

C  -  1  ^ 

x  "  jwe  3x3z 

and  since  there  is  no  z-variation,  E2X  vanishes.  Furthermore,  E^  '*‘,b,zi 

also  vanishes  for  z  <  0,  Hi  is  a  small  diffracted  field  for  z  <  0.  Hence, 

lx 


(28)  becomes 


«  a 


/  E'z(x',b,z)dz  =  -  /  {  /  E?7(x,b)H;y(x,b,z)dx}dz  (29) 


2z '  *  7  lx 


0  “  0  0 
Both  Ejz  and  H|x  have  the  same  z-variation,  but  E^z  has  no  z-variation.  Thus, 

(30) 


From  (23) 


Therefore, 


E‘z(x',b,z)  =  -/  E2z(x,b)H^(x,b,z)dx 


H'x(x,b,z)  =  y-  E[z(x,b,z) 
s 


From  (13b) 


E jz ( x ' ,b,z)  =  -/  y-  E;7(x,b,z)E?7(x,b)dx 


0  Zs  lz 


2z' 


-Y  Z 

1  m 


(31) 


Eu(x'y’z)  ”  ky  y™  s1n(kyy)dky 


jus0  s’"l  a  '  aySz  *ri/y) 
Hence,  (31)  becomes 


(32) 


i  5i„(22l)  3 
Jwe0  m=l  * 


a  '  3y3z  rm 


*Jy) 


y=b 


-o/Ss4  X  ^  **  *”<y)U  E2z(X,b’dX  <33> 
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Equating  term-by-term  for  each  m, 


sin( 


mux 

a 


(TVnXv  p 

a  'b2z 


(x ,b)dx 


(34) 


J 

- 1 ///rlrrtrtrr'  - 

t 

cr=oo  crm=oo  cr  =  co 


j 

>////?  ■>? T?  A?-/ / T7 rT7 7777 

CTm  =  0O 


(a) 


(b) 


Fig.  17.  (a)  Geometry  for  the  unperturbed  problem 

(b)  Approximation  of  (a). 


As  a  first  aproximation  to  find  E22,  the  fields  on  the  magnetic  strip  laid  on 
an  electrically  conducting  plane  is  approximated  by  the  fields  on  a  perfect 
magnetic  conductor  plane  (Fig.  17).  Then  E2z(x,b)  of  the  source  at  (x',b)  can 
be  found  using  the  potential  function 

*h  =-jjj  Hj2)(k|x-x'|)  (35) 
and 


H 


V 


X 


For  TMZ  fields. 


r  =  J-  (i-,  +  k2)^ 

Z  jue  '  3zZ 


SO 


E2z<x-b>  *  --h 


(36) 


Thus,  (34)  becomes 


m  th  mode  can  be  defined 


Assume  that  the  m=l  mode  is  predominant  in  the  open-walled  waveguide.  Then, 


2$  =  and  assume  that 


-Y1Z 

^(y)  =  A  cos(kyly)  e 


where  is  complex. 
Then , 


(42) 


Finally,  (39)  becomes 


—  k  ,  tan(k  .b) 
s  Jue0  yl  v  yl  ; 


(43) 


kyl  tan(kyib) 


a  a  iH 

/  /  e  a  (k|x-x' | )dx  dx' 

0  0  u 


(44) 


This  final  equation  has  a  Hankel  function  with  complex  argument  in  its 
integrand,  so  numerical  methods  must  be  used  to  solve  it  for  k^.  Then  Yj  is 

Y1  /  V  Kyl  K0 
and,  by  the  assumptions  previously  taken 


Solving  (44)  for  ep  from  the  measured  y  is  very  difficult  and  takes  much 
computing  time.  However,  computing  y  from  assumed  values  of  er  is  relatively 
simple.  It  is  then  possible  to  plot  a  map  of  er  as  a  function  of  y  as  shown 
In  Fig.  18.  Some  of  the  measured  y's  shown  in  Fig.  11  are  also  shown. 
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IV.  DISCUSSION 

The  three  data  points  which  are  located  below  the  =  0  curve  in  Fig.  18 
do  not  fit  well,  since  is  negative  (as  er  =  -  je^,  both  t'r  and  e" 

should  be  positive).  The  probable  reason  for  the  discord  of  the  numerical  and 
experimental  results  is  that  the  approximations  taken  in  the  theory  are  too 
rough,  especially  in  the  assumed  field  in  the  unperturbed  state,  (36).  For  a 
better  approximation ,  the  original  situation  of  Fig.  17(a)  can  be  approximated 
by  removing  the  magnetic  strip  and  replacing  the  source  with  a  source  of  2J 
(original  source  plus  its  image)  situated  between  two  thin  electrically 
conducting  half-planes  as  shown  in  Fig.  19. 


\ 

<J=00 


cr=co 


Fig.  19.  Approximated  problem  for  Fig.  17(a). 

In  this  case  the  potential  function  of  the  form 

*h=7T  l  *m  s1n(£f*)H£Z)(Mx-x,l)  («) 

J  m=l 

should  be  tried  instead  of  (35).  This  function  satisfies  the  boundary 
conditions  at  the  both  edges,  x=0  and  x=a,  of  the  conducting  half  planes. 
Again,  the  m=l  mode  would  be  dominant.  Hence, 

r-2jSln(^)H<2)(k|x-x'|) 


(46) 


So  instead  of  (36),  a  better  approximat i on  to  the  unperturbed  field  is 

2 

E2z(X’b)=“2^  sin(^)  H(2)(k|x-x'|)  (47) 

And  finally,  (44)  would  be  replaced  by 

k2 

kyl  tan (kylb)  =  j  — ^  /  /  sin 2(J1)h[)2^  (k  | x-x '  |  )dx  dx '  (48) 

which  should  improve  the  results  significantly. 

Another  problem  is  that  the  computed  y's  are  concentrated  in  a  narrow 
strip  for  >  4.0.  This  characteristic  makes  it  difficult  to  accurately 
find  e  from  the  measured  y.  even  if  the  theory  is  correct.  We  believe  that 
this  difficulty  can  be  overcome  by  filling  the  waveguide  with  a  low-loss 
dielectric  material  with  e'r  =  4  to  6.  Introducing  such  a  dielectric  in  the 
system  will  scale  down  the  ep  of  the  material  relative  to  that  inside  the 
waveguide.  This  would  make  the  most  of  y‘s  fall  in  the  wide-spread  region  of 
the  map.  Conversion  to  the  system  to  use  a  dielectric-filled  waveguide  is 
planned. 


V.  FUTURE  STUDIES 


The  accuracy  in  measuring  a  can  be  substantially  increased  by  using  a 
"phase-insensitive  homodyne  detector"  used  by  King  and  Yen  [1981]  which  is 
similar  to  a  quadrature  SSB  modulator  with  output  of  the  form 

K,  Iblm  cos(u)  t+$) 
h 1  '  '  m  ' 

instead  of 

| b  jm  cos($)cos(wmt) 

as  in  (2).  Hence,  a  continuous  variation  of  |b|  can  be  measured  to  determine 
a  in  the  phase-insensitive  system,  while  only  several  peaks  of  the  output  are 
used  in  the  conventional  phase-sensitive  system  as  seen  in  Fig.  10. 


Fig.  20.  Phase-insensitive  homodyne  system. 
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Moreover,  this  phase-insensitive  system  detects  one-way  transmission  of  the 
signal 

b  -  bQ  e-“z  e-jBZ 

rather  than  the  round-trip  backscattered  signal.  Thus,  this  system  is  twice 
as  sensitive  (as  measured  in  dB).  In  this  system  the  phase  gz  is  measured  by 
a  phase  meter. 


Fig.  21.  Comparison  of  the  outputs  of  phase-insensitive 
system  and  phase-sensitive  systems. 


Both  the  invasive  and  the  non-invasive  methods  used  have  some 
difficulties.  However,  most  of  them  can  be  reduced  using  different  devices 
for  coupling  the  fields  to  the  medium. 

For  the  Invasive  method,  the  surface  wave  launching  arrangement  shown  in 
Fig,  22  has  several  advantages  over  that  shown  in  Fig.  2. 
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Fig.  22.  New  invasive  surface  wave  structure. 

The  microwave  instrumentation  is  rigidly  attached  to  a  coaxial  line  which 
feeds  a  surface  wave  structure  through  a  tapered  transition.  The  surface  wave 
structure  consists  of  a  dielectric  sheath  on  the  stationary  metal  tube.  The 
sheath  is  terminated  with  a  gradual  taper  just  above  the  medium  interface,  and 
the  metal  tube  intrudes  to  a  depth  of  about  2A  into  the  medium.  A  modulated 
monopole  scatterer  mounted  at  the  end  of  a  plastic  tube  is  inserted  into  the 
metal  tube,  extending  into  the  medium  through  a  slot  cut  in  the  metal  tube. 
Modulation  of  the  PIN  diode  terminating  the  scatterer  Is  accomplished  by  leads 
running  through  the  plastic  tube  to  an  external  AF  oscillator. 

Advantages  of  this  new  system  are  that: 

a)  The  surface  wave  structure  tends  to  focus  the  radiation  into  the  medium 
along  the  metal  tube  so  that  the  field  intensity  is  greater  than  for  the 
old  system  In  Fig.  2. 


‘■H 
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b)  The  microwave  source  is  a  rigid  integral  part  of  the  probing  mechanism. 

c)  A  monopole  more  nearly  probes  the  field  at  a  point  (in  contrast  to  the 
area  of  a  loop).  It  is  easier  to  tune,  and  the  scattered  field  is 
generally  larger. 

d)  The  leads  being  used  to  modulate  the  PIN  diode  are  shielded  by  the  metal 
tube. 

e)  The  operating  frequency  range  is  much  broader  than  for  the  old  system. 
The  upper  frequency  limit  is  determined  oy  the  presence  of  higher  order 
modes  propagating  on  the  surface  waveguide,  and  the  lower  frequency  limit 
is  generally  determined  by  the  insertion  length  of  the  metal 
tube  (~  2x),  and  the  system  sensitivity. 

For  the  non-invasive  method  the  open-walled  rectangular  waveguide  has  a 
very  limited  operating  frequency  range.  Approximately,  this  range  corresponds 
to  cutoff  for  the  dominant  TEjq  mode  of  a  closed  waveguide  at  the  low  end,  and 
to  cutoff  of  the  next  higher  order  mode  at  the  high  end.  The  split  rigid 
coaxial  line  in  Fig.  23  has  a  cut  off  frequency  of  zero.  This  considerably 
extends  the  lower  end  of  the  frequency  range  until  the  physical  size 
limitation  is  reached  (again,  approximately  2X). 

Another  suggested  non-invasive  method  uses  coplanar  strip  lines  laid 
directly  on  the  medium  to  measure  the  propagation  constant  of  a  quasi-TEM 
wave.  Its  configuration  is  shown  in  Fig.  24.  Since  there  are  two  conductors 
and  two  dielectrics  (air  and  the  medium),  the  dominant  mode  is  quasi-TEM. 
Therefore,  its  cut  off  frequency  is  zero.  The  lower  end  of  the  usable 
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Fig.  23.  Split  coaxial  line  traveling  wave  sensing  head. 


bandwidth  will  chiefly  be  determined  by  the  length  of  the  strip  lines.  The 
higher  end  will  be  determined  by  appearance  of  any  higher  order  modes.  Hence 
this  method  gives  a  wide  operating  frequency  range. 

Another  advantage  of  this  method  is  that  it  can  adapt  to  a  slowly  varying 
uneven  surface.  The  coupling  of  the  wave  to  the  medium  can  be  increased  by 
increasing  the  width  of  the  gap  between  the  two  conductors.  Wide  conductors 
also  enhance  the  coupling  and  can  reduce  the  effect  of  surface  roughness  as 
suggested  by  Chang  and  Wait  [1974]. 

The  coplanar  strip  line  can  be  fed  using  a  semi-iigid  coaxial  line  by 
grounding  the  shield  to  one  strip  and  the  extended  center  conductor  to  other 
strip,  as  shown  in  Fig.  24.  The  line  can  be  terminated  in  ZQ  over  a  broad 
range  of  frequencies  using  an  electric  and/or  magnetic  absorbing  wedge. 

The  modulated  scattering  probe  can  be  either  an  electric  monopole  or 
magnetic  dipole,  terminated  in  a  PIN  diode  which  is  switched  at  The  probe 
can  be  moved  along  the  line  to  measure  a  and  6  of  the  traveling  wave. 
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Fig.  24.  Coplanar  strip  line  traveling  wave  sensing  head. 


For  a  half-buried  round  wire,  Chang  and  Wait  found  that 

1+e. 


T0"=  C~ 
J30 


r-jl/2 


which  means  that  the  effective  dielectric  constant  is  the  average  of  the  two 
half-spaces.  It  is  expected  to  t>at  the  measured  y  will  have  fairly  strong 
dependence  on  ep  of  the  medium  which  promises  more  accurate  measurement  of  er. 
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